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We experimentally tested the impact of browsing, defecation and urination corre- 
sponding to four different levels of moose population density on abundance and 
number of eco-morphological groups of soil mesofauna in a Swedish boreal fo- 
rest. The study was carried out in three fenced exclosures representing different 
levels of productivity (from nutrient-rich to nutrient-poor). Moose impact on soil 
communities was evaluated by analysing abundance, richness and diversity of 
mesofaunal groups, and by means ofa multi-taxa index based on ecomorphologi- 
cal groups, the Biological Soil Quality index (QBS). There was a negative impact 
of high moose densities on the abundance and richness of soil mesofauna. Fur- 
thermore, low-moderate moose densities had a positive impact on abundance and 
number of eco-morphological groups. We conclude that moose can have a large 
impact on soil-living mesofauna, and therefore probably soil decomposition pro- 
cesses and biological soil quality, in young boreal forest stands. 
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1. Introduction 


The selective feeding of large mammalian herbi- 
vores can significantly change plant community 
structure, species composition and productivity 
(e.g. Pastor & Naiman 1992, Hobbs 1996, Au- 


gustine & McNaughton 1998, Danell et al. 2006). 
These changes can further indirectly affect the 
abundance and diversity of animal groups at dif- 
ferent trophic levels and fundamental ecosystem 
processes such as the quantity and quality of 
litterfall and soil decomposition and mineraliza- 
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tion (Pastor & Naiman 1992, Augustine & Mc- 
Naughton 1998, Cote et al. 2004, Danell et al. 
2006, Suominen et al. 2008). 

The moose (Alces alces Linnaeus) has been 
shown to have a significant impact on the ecosys- 
tem of boreal forests, especially at higher popula- 
tion densities (e.g. Pastor & Naiman 1992, Danell 
et al. 1994, Persson et al. 2007, Suominen et al. 
2008). The population densities in some areas of 
Sweden are among the highest in the world and 
constitute densities probably not experienced in 
post-glacial times (Cederlund & Bergström 
1996). Moose mainly affect the ecosystem by se- 
lective browsing in the tree and shrub-layers in 
young forest stands. Young trees and shrubs make 
up a major part of their diet both in summer and 
winter (Cederlund et al. 1980, Bergström & Hjel- 
jord 1987), and moose browsing at high popula- 
tion densities generally results in decreased tree 
height, a more open tree canopy, and increased 
dominance of less preferred species (McInnes et 
al. 1992, Pastor & Naiman 1992, Persson et al. 
2007). 

The changes in the tree and shrub layers have 
further been shown to result in lower litter pro- 
duction and subsequent lower soil nutrient levels 
at higher moose densities (Pastor & Naiman 
1992, Persson et al. 2005a), which can result in 
lower decomposition rates and soil nutrient cy- 
cling (Pastor & Naiman 1992, Pastor et al. 1993). 
Browsed trees can also be expected to reduce al- 
location of resources to root growth (Côte et al. 
2004) and thus affect mycorrhizal fungi. How- 
ever, such impacts on the ecosystem also strongly 
depend on moose population density, and moose 
at lower densities can even have positive effects 
on the above processes (Pastor & Naiman 1992, 
Persson et al. 2005a,b, 2007). For example, 
moose also deposit faeces and urine, which con- 
tain easily available soil nutrients which may 
have important effects on vegetation and soil mi- 
crobial biomass (as demonstrated experimentally 
in reindeer Rangifer tarandus; Van der Wal et al. 
2004), but the effects of selective feeding have 
been shown to be more important at larger spatial 
and temporal scales (Pastor & Naiman 1992). 

The soil fauna is important in maintaining 
ecosystem functioning, and fundamental soil pro- 
cesses may be maintained or enhanced by con- 
serving the biological diversity of soil ecosys- 
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tems (Stork & Eggleton 1992, Filser et al. 1995, 
Bengtsson et al. 2000). Furthermore, much of the 
biodiversity of forest ecosystems resides in soil 
(Behan-Pelletier & Newton 1999). Soil meso- 
fauna constitute an important animal group in the 
boreal forest in terms of ecosystem processes 
such as decomposition of organic matter, nutrient 
mineralisation and plant growth (Liiri et al. 
2002). They are linked with catalysis of organic 
matter decomposition and nutrient cycling, they 
have an important regulatory role in soil micro- 
bial processes and community structure, soil 
structure and hydrologic fluxes, and they acceler- 
ate plant litter decomposition through interac- 
tions with the microflora (Moore & Walter 1988, 
Wardle et al. 1999). 

Soil quality has been defined as “the function- 
ality ofa soil in its own environment, the capabil- 
ity to sustain plants and animal productivity, and 
to maintain or improve the air and water quality” 
(Karlen et al. 1997). Certain groups of soil inver- 
tebrates are sensitive to variations in the terres- 
trial environment, such as nutrient levels and pol- 
lution, and thus are good indicators of changes in 
soil quality (Van Straalen 1997, Ponge et al. 
2003). Therefore, soil zoocoenoses have been re- 
cognised as important descriptors of environmen- 
tal quality (Parisi 2001, Parisi et al. 2005), and 
some soil invertebrates in particular are consid- 
ered important bioindicators of soil quality. The 
results from research on the direct and indirect in- 
teractions between moose, the plant community 
and soil decomposition and nutrient cycling (e.g. 
Pastor et al. 1993, Persson et al. 2007, Persson et 
al. 2009) indicate that moose have a strong poten- 
tial to affect abundance and diversity of soil-liv- 
ing organisms. However, soil mesofaunal sensi- 
tivity to impacts of moose, or any other native fo- 
rest-living large herbivore, has as far as we know 
never been examined, although research on im- 
pacts of non-native herbivores has shown some 
significant impacts on soil mesofauna in addition 
to other measures of soil quality (Wardle et al. 
2001). 

In order to examine moose impacts on soil 
mesofauna, we experimentally tested how simu- 
lated browsing, defecation and urination corre- 
sponding to four different population densities of 
moose (including the scenario with no moose 
present) affected various measures of soil meso- 
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fauna abundance, richness and diversity in sites 
ranging from very nutrient-poor to nutrient-rich 
in a Swedish boreal forest. We used an experi- 
mental set-up which has been extensively studied 
to examine moose impact on the boreal ecosys- 
tem, particularly on primary productivity and 
litterfall (Persson et al. 2005a,b, 2007). Based on 
previous findings from the experiment (as sum- 
marised above), we predicted that high levels of 
simulated moose density would have a negative 
impact on soil mesofauna, whereas low-moderate 
moose densities would have a positive impact. 
We further tested if abundance, richness and di- 
versity of soil mesofauna groups were correlated 
with litter production. Litter is an important factor 
regulating soil mesofauna (Blair et al. 1988, 
Eaton et al. 2004) and quality and quantity of de- 
tritus have been shown to affect mesofauna densi- 
ties (Seasted 1984, Takeda 1987). The decrease in 
litterfall with increasing moose densities previ- 
ously documented from the experimental set-up 
(Persson ef al. 2005a) suggests that abundance 
and diversity of soil mesofauna should be posi- 
tively correlated with litter production. 


2. Materials and methods 
2.1. Study sites 


The study was performed in an established exper- 
imental field set-up in the middle boreal zone of 
northern coastal Sweden situated 50-90 km N 
and NW of Umea (63°50’N, 20°18’E), with 
fenced 70 x 70 m exclosures (Fig. 1) situated in 
young forests dominated by Scots pine (Pinus 
sylvestris L.) and birch (Betula pubescens Ehrh. 
and B. pendula Roth.) with more sparsely occur- 
ring aspen (Populus tremula L.), rowan (Sorbus 
aucuparia L.) and willows (Salix spp.). The sites 
were chosen to represent the full range of the hab- 
itat productivity gradient in the study area: 1) a 
nutrient-poor lichen-rich pine heath (Légda- 
berget); 2) a site dominated by bilberry (Vaccini- 
um myrtillus L.) and lingonberry (V. vitis-idaea 
L.) in the shrub layer (Ralberget), and a mixture 
of birch, pine and Norway spruce (Picea abies 
(L.)) in the tree layer; and, 3) a very nutrient-rich 
site with a high intermixture of deciduous trees, 
ferns and forbs (Mörtsjöstavaren). For conve- 
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Fig. 1. The experimental design with fenced 
exclosures (4,900 m’), one exclosure per each of 
three habitat types (i.e. sites). Inside fenced 
exclosures we simulated browsing, defecation and uri- 
nation corresponding to 3 different moose densities in 
treatment plots of 25 x 25 m. One plot was left un- 
touched to represent no moose, and constituted the 
statistical control plot. There were buffer zones of 5 m 
between treatment plots as well as between treatment 
plots and the fence. Low-moderate moose density 
corresponds to 1 moose per km’ in winter, high to 3 
per km’ and extremely high to 5 per km’. 


nience, we use the term ‘habitat type’ to refer to 
the productivity gradient. 


2.2. Experimental design 


Full details of the experimental set-up, and the 
protocols to simulate browsing, defacation and 
urination have been published elsewhere (Pers- 
son 2003, Persson et al. 2005a), so we provide 
only a summary here. 

In each exclosure, browsing, defecation and 
urination corresponding to 4 levels of moose po- 
pulation density, from no moose present (the sta- 
tistical control) to one of the highest local densi- 
ties experienced in Sweden (5 moose per km’ in 
winter density) had been simulated since 1999 in 
plots of 25 x 25 m within each exclosure (Fig. 1). 
Protocols to simulate moose browsing were de- 
veloped to match as closely as possible known 
behaviour, diet and intake rates of moose in the 
area (Persson et al. 2000, Persson et al. 2005a, 
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Table 1. — a. Soil mesofaunal groups (based on Parisi et al. 2001) used in the analysis. — b. Abundance (individ- 
uals per 3,000 cm’ soil volume unit) for each relevant eco-morphological group identified, shown per each 
exclosure subject to different simulated moose densities and nutrient levels. — c. Total abundances (from part b) 
as well as diversity indices and the Biological Soil Quality index (QBS), shown per each exclosure subject to dif- 
ferent simulated moose densities and nutrient levels. 


1 





(a) Eco-morphological groups 





Group name 


Characteristics (where relevant) 





Collembola1 
Collembola 2 
Collembola 3 
Collembola 4 
Collembola 5 
Collembola 6 
Collembola 7 
Cicada 
OHemip 
Thysanoptera 
Coleoptera’ 
Coleoptera2 
Formicidae 
OHymen 
DipL 

Hol-L 

Hol-A 
Protura 
Diplura 
Opiliones 
Aran-L 
Aran-S 
Acarina 
Chil-L 

Chil-S 


Epigeous forms, complex pigmentation, well developed appendages 
Epigeous forms not associated with grass, trees, shrubs 

Small, modest pigmentation, average appendage development, 
Hemi-edaphic, not elongated appendages, cuticle pigmented 
Hemi-edaphic, reduced ommatidia, scarcely developed appendages 
Eu-edaphic, unpigmentated, furca reduced, ommatidia reduced 
Eu-edaphic, unpigmentated, furca absent, short appendages 
Larvae 

Other Hemiptera (epigeous or root-feeding forms) 


Adults, epigeous, >2mm 
Adults, epigeous, <2mm 


Other Hymenoptera 

Diptera larvae 

Other holometabolous insects, larvae 
Other holometabolous insects, adults 


Araneae >5mm 
Araneae, <5mm, scarcely pigmented 


Chilopoda >5mm, well developed legs 
Chilopoda <5mm 





(b) Soil mesofauna abundance 





Moose dens.# 
Productivity $ 
Collembola1 
Collembola2 
Collembola3 
Collembola4 
Collembola5 
Collembola6 
Collembola7 
Cicada 
Ohemip 
Thysanoptera 
Coleoptera’ 
Coleoptera2 
Formicidae 
Ohymen 

DipL 

Hol-L 

Hol-A 

Protura 
Diplura 


jo) 


E EE E fo 
OWNND =~ NCH 


oooff HAH WOCO-O0 


1 3 5 0 í 3 5 0 1 3 5 
E L L M M M M H H H H 
0 24 11 0 6 12 13 21 1 12 7 
8 28 4 12 15 0 4 2 38 24 0 
23 6 4 85 23 57 21 82 61 11 13 
24 22 5 4 1 3 25 25 45 14 24 
58 6 11 37 8 1 15 24 41 44 48 
7 26 7 39 156 22 8 15 46 91 74 
9 15 0 15 33 18 4 89 211 49 14 
1 0 0 0 1 0 0 2 3 0 0 
1 1 1 2 9 3 5 5 5 3 0 
0 0 0 0 1 1 0 2 2 0 0 
1 1 0 0 2 0 0 0 0 0 2 
0 0 0 3 0 0 3 1 2 0 0 
0 0 0 2 13 0 0 7 1 1 5 
1 0 0 0 0 0 0 1 0 0 0 
3 8 1 7 12 1 1 0 18 12 19 
1 19 6 13 26 55 4 21 16 26 7 
3 1 0 6 1 0 1 0 0 1 1 
0 0 0 3 2 1 0 0 1 1 0 
0 0 0 0 1 0 0 1 3 2 0 
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Table 1, continued 





Elia et al. + ENTOMOL. FENNICA Vol. 21 








Opiliones 0 0 0 0 3 A 0 0 1 0 1 0 
Aran-L 0 0 0 0 0 0 0 0 1 1 0 0 
Aran-S 0 0 2 0 0 1 1 3 1 4 0 0 
Acarina 623 797 419 349 734 112 384 164 928 687 728 512 
Chil-L 1 0 0 0 0 0 1 0 1 1 1 0 
Chil-S 1 1 0 0 0 0 0 0 0 0 0 0 
(c) Total abundance, total taxa and derived indices 

Moose dens.# 0 1 3 5 0 1 3 5 0 1 3 5 
Productivity $ L L L L M M M M H H H H 
Tot. abundance 704 936 578 399 965 424 560 271 1,230 1,187 1,021 726 
Total taxa 15 15 14 10 15 20 14 14 20 20 17 12 
Menhinick 0.57 049 058 0.50 048 0.97 059 O85 0.57 0.58 053 0.45 
Margalef 5.27 5.05 5.07 3.84 503 7.61 509 575 647 651 565 4.19 
Simpson 0.78 0.73 0.54 0.77 0.59 0.22 049 039 0.58 0.37 0.52 0.51 
Shannon 0.61 0.69 1.19 063 1.02 1.94 1.16 153 1.06 1.50 1.19 0.41 
QBS 97 94 68 51 102 134 97 82 119 138 127 67 





# Individuals per km? 
$ L = low, M = moderate and H = high 


2007), which varied according to plant species 
and season. The removal of biomass to simulate 
moose browsing in the treatment plots was calcu- 
lated for the different moose densities based on an 
estimated daily food intake of moose of 5 kg dry 
mass in winter and 10 kg in summer (Persson et 
al. 2000), and a winter and summer season of 180 
days each. The species composition of food 
plants in the harvest was based on the composi- 
tion of moose diet in boreal habitats in Sweden 
(Cederlund et al. 1980, Bergström & Hjeljord 
1987) and included trees (pine, birch, aspen, wil- 
low and rowan) and shrubs (bilberry, heather 
Calluna vulgaris (L.) Hull and fireweed Epilo- 
bium angustifolium L.). 

Browsing was simulated by clipping living 
shoots from a range of known food plants at ap- 
propriate heights and in appropriate dietary pro- 
portions. The height at which vegetation was 
sampled was also varied according to season 
(vegetation was usually under snow up to a height 
of 50 cm in winter). The mass of all removed 
plant material was measured each time to 0.1 g 
precision. Clipping was carried out until the esti- 
mated biomass for the moose density in question 
was obtained. All removed plant material was de- 
posited outside the study area. In the growing sea- 
son (May through October) clippings were taken 
once a month. For winter, when plants are dor- 


mant, clippings for the whole winter were taken 
in April before the growing season started. We 
also added natural moose dung (collected at a 
nearby moose farm where moose were fed a natu- 
ral diet; Nyberg & Persson 2002) and urine (arti- 
ficially made by dissolving urea in water, the con- 
centration of which was based on moose studies 
in North America; Persson 2003) in proportion to 
the simulated moose density. 


2.3. Soil sampling 


In 2004, after 5 years of experimental treatment, 
soil samples were collected. In the mid-growing 
season (July), 3 soil cores per treatment were ran- 
domly collected within each exclosure. Sampling 
was carried out by driving iron cores into the soil 
in dry weather. Each soil core was square in 
cross-section 10 cm deep and 100 cm’ in area; 
thus 1,000 cm’ soil was collected per sample, and 
then used for extracting soil arthropods via modi- 
fied Berlese-Tullgren funnels. Data on litter fall 
were available from a previous study (Persson et 
al. 2005a), where we collected plant litter from 
the tree and shrub layers for one year (i.e. one lit- 
ter cycle). Thirty traps (plastic pots of height 35 
cm, top diameter 50 cm and bottom diameter 35 
cm) were placed in each treatment plot in all 


ENTOMOL. FENNICA Vol. 21 * Impact of moose on soil fauna 


173 


Table 2. Test statistics for the ANOVA models of the effects of treatment (i.e. level of simulated moose density) 
and habitat type (block factor) on total no. of mesofauna individuals (Tot. mesofauna), total no. of eco-morpholog- 
ical groups (Tot. groups), total no. of Collembola and total no. of Acari sampled per treatment area on a 3,000 
cm’ soil volume unit, as well as the Biological Soil Quality index (QBS). F-values, degrees of freedom and p-val- 
ues are presented for the full two-way ANOVA, and for treatment and site effects separately. 








Variable tested Statistic Full model Treatment Habitat type 
df= 5,6 df = 3,6 df= 2,6 

Tot. mesofauna F 7.08 5.20 9.91 

p 0.02* 0.04* 0.01* 
Tot. groups F 6.00 6.97 4.54 

p 0.02* 0.02* 0.06 
QBS F 9.17 9.85 8.15 

p 0.01* 0.01* 0.02* 
Tot. Collembola F 9.04 5.51 14.35 

p 0.009* 0.04* 0.005* 
Tot. Acarina F 3.51 2.93 4.38 

p 0.08 0.12 0.07 





* p values reporting statistical significance (p < 0.05) according to ANOVA. 


exclosures, for a total of 120 per exclosure. The 
traps were emptied monthly during the growing 
season and as soon as possible after snowmelt in 
spring. The litter was dried at 40° C, and produc- 
tion of litter (g dry mass per m? per year) for each 
exclosure and treatment plot was calculated. 


2.4. Soil mesofauna descriptors 


Arthropods were counted by means of micro- 
scopic observations and classified according to 
the biological form approach (Sacchi & Testard 
1971). Biological forms were separated into 
groups having homogeneous morphological 
characters (listed in Table 1), according to the 
classification of Parisi et al. (2005). Total meso- 
fauna abundance, number of eco-morphological 
groups and total number of Collembola and 
Acarina were counted by pooling together the 
three samples collected in every treatment area. 
In addition, we consider an indicator of bio- 
logical soil quality, the Biological Soil Quality in- 
dex (QBS; Parisi 2001, Parisi et al. 2005). The 
QBS is intended to assess soil conditions via anal- 
ysis of microarthropod communities and is based 
on the concept that a higher soil quality is re- 
flected by a higher number of microarthropod 
groups well adapted to the soil habitat (Parisi et 
al. 2005). The QBS takes into account micro- 
arthropod adaptation to the soil by giving a nu- 


meric value to each biological form (sensu Sacchi 
& Testard 1971) in relation to its degree of adap- 
tation to the edaphic environment, e.g. reduction 
or loss of pigmentation and visual apparatus, 
streamlined body form and reduced appendages. 
Higher scores indicate a greater level of adapta- 
tion to soil life for a given mesofaunal group. For 
example, groups such as Symphyla, Acarina and 
Protura have a value of 20, Diptera larvae 10. 
Other taxa, such as Collembola, encompass a 
number of groups with a range of scores (from 1 
to 20) reflecting their diversity of adaptation to 
the edaphic environment (further details given in 
Parisi et al. 2005). The QBS index is the sum of 
the scores across the defined groups. 

Given the importance of the biodiversity of 
soil systems to fundamental soil processes (Stork 
& Eggleton 1992, Filser et al. 1995, Bengtsson et 
al. 2000), we also expect impacts of simulated 
moose density on diversity measures. The eco- 
morphological groups classified by the QBS in- 
dex were therefore also used to calculate indices 
describing different aspects of soil fauna diversity 
(Magurran 2004): species richness — Menhinick 
(Menhinick 1964) and Margalef (Margalef 
1957); dominance — Simpson (Simpson 1949); 
and, diversity — Shannon-Weaver (Shannon & 
Weaver 1948). All indices were calculated as a 
single value for each replication of treatment (i.e. 
by pooling together counts from the three treat- 
ment cores collected in every treatment area). 
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2.5. Statistics 


Effects of different levels of simulated moose po- 
pulation density on soil mesofauna abundance 
(number of individuals per treatment), group 
richness (number of groups per treatment, as de- 
fined in Table 1), diversity (Menhinick, Margalef, 
Simpson and Shannon-Weaver) and QBS index 
were analyzed using analysis of variance 
(ANOVA) models with level of simulated moose 
density as the nominal experimental factor (NV = 
4) and site (i.e, habitat type) as a block factor. Ef- 
fects on abundance were first analyzed for all soil 
mesofauna pooled and thereafter for the domi- 
nant groups Collembola and Acarina separately. 
When the ANOVA was significant (p < 0.05), 
post-hoc Tukey tests were carried out to show 
which specific treatments differed from each 
other. Relationships between litter production 
and mesofauna abundance, richness, diversity 
and QBS index were analysed using linear re- 
gression. 


3. Results 


A total of 9,001 individuals of soil invertebrates 
were extracted and classified according to the 
groups in Table 1. The ANOVA models showed 
that total number of organisms (i.e. total abun- 
dance) and the total number of eco-morphologi- 
cal groups sampled per treatment area on a 3,000 
cm’ soil volume, as well the QBS index, were sig- 
nificantly affected by both simulated moose den- 
sity (treatment) and habitat type (Table 2). 

The total number of groups per treatment and 
the QBS index showed the same response pattern: 
the highest value in the plots where we had simu- 
lated low-moderate moose density (1 per km’ in 
winter density), and the lowest value in the plots 
where we had simulated a very high moose den- 
sity (5 moose per km’) (Fig. 2). For total abun- 
dance of soil mesofauna, the response was more 
linear, with a steady decrease along the gradient 
and significant differences between no simulated 
moose density and very high simulated density 
(Fig. 2b). Of the two numerically most important 
mesofauna groups, Collembola and Acarina re- 
spectively (Table 1), the total number of Collem- 
bola sampled per treatment area was highest at 
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Fig. 2. The effect of different levels of simulated 
moose population density. — a. Total no. of eco-mor- 
phological groups, averaged over sites (n = 3). — b. 
Total no. of soil-living arthropods, averaged over sites. 
—c. Biological Soil Quality index (QBS). — d. Total no. 
of Collembola. Only variables showing significant (p < 
0.05) effects of simulated moose density are shown. 
Means and standard errors are presented. Different 
letters above the columns indicate that treatments are 
significantly different from each other at the p < 0.05 
level (Tukey test). Estimates of total mesofauna abun- 
dance and diversity are per volume soil sample (3,000 
cm). 


low/moderate moose density and lowest in the 
very high density treatment (Fig. 2d). For total 
number of Acarina, there was no significant ef- 
fect of treatment. However, both total number of 
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Fig. 3. The effect of habitat type, corresponding to dif- 
ferent nutrient levels. — a. Total no. of soil-living arthro- 
pods. — b. Simpson index. — c. Biological Soil Quality 
index (QBS). — d. Total no. of Collembola. Means av- 
eraged over treatments (n = 4). L = low nutrient level, 
M = intermediate nutrient level, H = high nutrient level. 
Other details as in Fig. 2. 


Collembola and Acarina showed a similar gen- 
eral pattern to the other measures, with the lowest 
number in the plots for which we had simulated a 
very high moose density (Fig. 2d, Table 1). 
ANOVA tests for diversity measured by the other 
biotic indexes found no significant effect of treat- 
ment on any index (F, , < 2.39, p > 0.17). 

Site effects were significant for total meso- 
fauna abundance, total number of Collembola, 


the QBS (Table 2) and the Simpson index (F). 


Values were higher in the high nutrient sites with 
the exception of the Simpson index, where sites 
with intermediate nutrient levels had the lowest 
diversity (Fig. 3). Differences approached signif- 
icance (p = 0.06—0.07) for total number of groups 
and total number of Acarina (Table 2), in each 
case the highest value being in the high nutrient 
site. 

As shown in Fig. 4, there were positive linear 
relationships between litter production (g dry 
mass per m° per year) and total number of soil me- 
sofauna (F, ,, = 16.70, p = 0.003, r = 0.63), total 
number of eco-morphological groups (F,,, = 
7.43, p = 0.02, 1 = 0.43), the QBS index Fo = 
6.82, p = 0.026, r = 0.41), total number of 
Collembola (F, ,, = 29.79, p < 0.001, r = 0.75) 
and total number of Acarina (F,,, = 4.51, p = 
0.059, 7? = 0.31). 
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4. Discussion 


Our results show that moose can indirectly have a 
significant impact on the abundance and number 
of eco-morphological groups of soil mesofauna, 
and on QBS index, in young boreal forest stands, 
and that the effect is strongly dependent on moose 
population density. The main conclusion is that, 
despite the wide range in habitat productivity, the 
effect of simulated moose browsing was consis- 
tent across the productivity gradient. There was 
also a general positive effect of nutrient level on 
mesofauna abundance and QBS as expected, al- 
though Simpson’s index showed the lowest val- 
ues at intermediate nutrient levels, possibly sug- 
gesting non-linear responses of certain diversity 
measures to productivity gradients. 

The number of groups, the abundance of 
Collembola and QBS index all showed the high- 
est value at low-moderate levels of moose popu- 
lation density indicating that moose at normal 
densities (around 1.0 moose per km’ for many ar- 
eas in Sweden in winter) can actually have a posi- 
tive effect on soil mesofauna in young boreal fo- 
rest stands. However, at higher moose densities, 
corresponding to about 3 moose per km? or more, 
moose can have a negative impact on soil arthro- 
pods. Densities of 2 — 3 moose per km’ have oc- 
curred over large areas in Sweden (Cederlund & 
Bergström 1996) and the local density (e.g. in 
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Fig. 4. Regressions between litter production (g dry 
mass per m per year) and soil variables. — a. Total 
number of eco-morphological groups. — b. Total meso- 
fauna abundance. — c. Biological Soil Quality index 
(QBS). — d. Total no. of Collembola. — e. Total no. of 
Acarina. Results of statistical tests are given in the 
text. 


winter concentration areas) can be considerably 
higher (up to 5 moose per km? or more). Although 
the general patterns in mean abundance in 
Acarina were similar to those for the number of 
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groups, the abundance of Collembola and QBS 
index, suggesting some association with lower 
moose density sites (Table 1), treatment effects 
were not significant. This is probably due to the 
fact that we considered the abundance of all soil- 
dwelling Acarina sampled per treatment area (see 
below). 

The QBS was formulated as an indicator of 
biological soil quality (Parisi et al. 2005). Whilst 
previous work has shown the QBS to be sensitive 
to broad land use and farming system, so far there 
is little evidence to suggest that the QBS is closely 
related to broader measures of soil quality such as 
pollution (Parisi et al. 2005) or organic carbon 
content (Gardi et al. 2002a). In this study, the 
QBS index proved to be a sensitive tool to esti- 
mate the simulated effects of moose density on 
soil mesofauna, showing significant differences 
when the other main diversity indices (Shannon- 
Weaver, Menhinhick, Margalef and Simpson) 
failed to detect significant effects. In addition, the 
QBS showed significant associations with nutri- 
ent levels and a positive relationship with litter 
quantity. As there are also known negative im- 
pacts of high moose densities on wider aspects of 
soil quality from the same experimental set-up, 
namely C and N contributions (Persson et al. 
2005a), we argue that these results suggest that 
the QBS may be a useful indicator of biological 
soil quality. Nevertheless, there may be some im- 
provements that could be made in the compo- 
nents of the index, in particular Acarina are con- 
sidered as a single group, thus pooling together a 
(probably high) number of species having high 
ecological diversity (Coleman et al. 2004), differ- 
ent feeding habits and using a wide range of eco- 
logical strategies (Vannier 1985, Lavelle & Spain 
2002). Further development of the QBS index 
should consider sub-division of some of the key 
groups, the Acarina in particular. 

It is difficult to draw firm conclusions on the 
more general impact of moose on soil mesofauna 
with only one sampling period and limited num- 
ber of samples. However, the fact that our results, 
despite the limited sampling and the large range 
in habitat productivity, showed significant effects 
on several variables strongly suggests that moose 
browsing is very likely to be an important factor 
determining pedofaunal communities in young 
boreal forest stands. The effects of moose at dif- 
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ferent population densities on the soil mesofauna 
are likely to be caused by an indirect effect medi- 
ated by a chain of steps. The most important un- 
derlying mechanism is probably the effect of 
moose browsing on the growth of trees and 
shrubs and consequently on canopy structure and 
litter production (McInnes et al. 1992, Pastor & 
Naiman 1992, Persson et al. 2005a, 2007). 

The highest values in soil mesofauna abun- 
dance (mainly a result of the response in the dom- 
inant group Collembola), number of groups and 
the QBS index at low-moderate moose density 
might be explained by a slightly more open can- 
opy increasing the amount of light reaching the 
forest floor and supporting more luxuriant herb- 
layer vegetation (McInnes et al. 1992, Persson et 
al. 2007). Previous studies (e.g. Atlegrim & Sjö- 
berg 2003) indicate that the abundance of terri- 
colous Collembola in boreal forests increases 
from shaded to more open forests. At the same 
time, the browsing pressure was probably too low 
to have any significant negative effects on annual 
litter production from the tree and shrub layers at 
low-moderate densities (Persson et al. 200Sa). At 
higher moose densities, it is likely that the large 
decrease in litter production (Persson et al. 
2005a) was a major explanation for the decrease 
in mesofauna abundance, number of groups and 
the QBS index, which is reflected in the signifi- 
cant positive correlations between litter and 
mesofaunal abundance and QBS (Fig. 4). Intense 
browsing substantially reduces both the amount 
of litter and nutrients returned to soil (Pastor & 
Naiman 1992, Pastor et al. 1993, Persson et al. 
2005a) which limits soil nutrient resources avail- 
able for ground-dwelling fauna. The linear re- 
gressions showing positive effects of litter fall 
thus provide a strong indication that the effects of 
the simulated moose browsing on litter produc- 
tion were a main factor mediating the observed 
effects on the soil mesofauna. This confirms that, 
with respect to high moose densities, the limita- 
tion in soil nutrient resources caused by browsing 
and therefore litter reduction may have a stronger 
impact on soil mesofauna than the return in dung 
and urine (Pastor et al. 1993). 

Collembola appeared to be more sensitive 
than other arthropods to moose impacts, showing 
a peak in abundance at moderate levels of moose 
density and low abundance at the highest densi- 
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ties. Collembola populations probably benefitted 
at first from a higher rate of sunlight reaching the 
ground and the associated increased decomposi- 
tion rate (Pastor et al. 1993, Augustine & Mc- 
Naughton 1998), thereafter decreasing when the 
impact of reduced litter production became 
greater than the positive effect of a more open 
canopy. Moreover, browsing by moose has been 
shown to have a significant negative effect on soil 
fungal communities and in particular on mycor- 
rhizal infections of fine roots of preferred browse 
species (Rossow etal. 1997, Persson et al. 2000). 

The fact that intense browsing substantially 
reduced the amount of leaf litter, likely also re- 
sulted in limited colonization by saprophytic 
fungi which are the most important decomposers 
in forest soils (Flanagan and Van Cleve 1983) and 
have been shown to be highly palatable to eda- 
phic Collembola known to be primarily fungi- 
vorous (Fitter & Sander 1992, Hopkin 1997, 
Rusek 1998, Jonas et al. 2007). 

At the same time, it is likely that high brows- 
ing pressure decreased allocation of resources to 
root growth (Côte et al. 2004), and thus also had a 
negative impact on mycorrhizal fungi, which is 
another food source for Collembola (Setälä 1995, 
Jonas et al. 2007). 

Since not only Collembola, but a range of soil 
invertebrates, feed on fungal spores and hyphae 
(Rabatin & Stinner 1988, Shaw 1992, Mc- 
Gonigle 1997), we suggest that an indirect nega- 
tive effect of intense browsing on soil fungi may 
have contributed to the pattern of mesofauna 
abundance by limiting their food sources. Thus, 
the negative impact on soil mesofauna at high 
moose densities might be due both to reduced lit- 
ter fall and reduced presence of saprophytic and 
mycorrhizal fungi, although more studies are 
needed to separate the effects of these variables. 

We are aware that moose can affect boreal 
ecosystems not only through their selective feed- 
ing, defecation and urination, but also through 
trampling, which has not been simulated in our 
experiment. Trampling certainly is a disturbance 
factor affecting composition of ground vegeta- 
tion, development of fine roots, mineralization 
rates of nitrogen, soil compaction (Persson et al. 
2000) and the movement of soil organisms 
through soil pore space (Wallwork 1970) and it is 
therefore likely to affect the diversity and density 
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of soil mesofauna in forest soils, as demonstrated 
in the study carried out by Battigelli et al. (2004). 
However, research has shown that browsing and 
thus biomass loss is the most significant impact of 
moose in boreal forest ecosystems, and, consider- 
ing the significant effects of decreased litter fall 
on several parameters, including vegetation rich- 
ness and nutrient availability (Suominen et al. 
2008, Persson et al., in press), it has been sug- 
gested that selective feeding affects the forest 
habitat more than trampling (Pastor & Naiman 
1992, Pastor et al. 1993, Suominen et al. 1999). 
Furthermore, the field work necessitated some 
(albeit light) trampling; the time spent in each 
treatment plot was in accordance with the differ- 
ent simulated moose densities. Further research 
on how moose trampling affects soil mesofauna 
in boreal forests is therefore needed, and will 
probably reveal complex effects and interactions 
with the other disturbance factors. 


5. Conclusions 


The importance of soil mesofauna for the integ- 
rity and functioning of terrestrial ecosystems is 
well recognised (Pankhurst 1997, Behan- 
Pelletier & Newton 1999). Especially important 
is the role that soil fauna dynamics have on soil 
organic matter degradation in boreal acid soils 
(Kurcheva 1960). We thus conclude that moose 
density can be a major factor affecting soil-living 
arthropods, and consequently early stage decom- 
position of soil organic matter and biological soil 
quality in young boreal forest stands. However, it 
is clear that further research is needed to fully un- 
derstand the processes involved. For moose man- 
agement, and for conservation of biodiversity and 
fundamental ecosystem processes, our results 
strongly suggest that impact of large herbivores 
on soil invertebrates should be taken into account. 
This is particularly important as moose popula- 
tions are increasing (Cederlund & Bergström 
1996) and densities in some areas are approach- 
ing those simulated in this experiment (5 individ- 
uals/km’). Such high moose densities are likely to 
be detrimental to soil quality, and hence manage- 
ment to maintain moose densities to low-moder- 
ate levels, where effects on soils are beneficial, 
should be encouraged. 
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